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Abstract. We investigate strong mechanical feedback for the single-electron 
tunneling device coupled to an underdamped harmonic oscillator in the high-frequency 
case, when the mechanical energy of the oscillator exceeds the tunnel rate, and for weak 
coupling. In the strong feedback regime, the mechanical oscillations oscillated by the 
telegraph signal from the SET in their turn modify the electric current. In contrast to 
the earlier results for the low frequencies, the current noise in not enhanced above the 
Poisson value. 
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1. Introduction 

One of the key questions in the field of nanoelectromechanical systems (NEMS) [T] 
is the effect of mechanical motion on the electric properties of the systems. It is 
equally important for undestanding of fundamental questions such as behavior of 
non-equilibrium, dissipative, driven systems and for the prospects of NEMS practical 
applications, such as switches, relays, and actuators. Electron motion in both external 
ac electromagnetic fields (see e.g. Ref. [2]) and subject to time-dependent noise [3] 
have been extensively studied in the past. The new feature brought by NEMS is that 
not only electrons move in the field created by the mechanical vibrations (for brevity, 
phonons), but also the phonons are created by the electron motion through the NEMS 
device. The phonons are generally driven out of equlibrium thus creating the feedback 
for the electron motion. 

In this Article, we concentrate on NEMS operating in the single electron tunneling 
(SET) regime. These include experimental realizations based on single molecules pE], 
semiconductor beams [5j EJ [7] , and suspended carbon nanotubes [H [H [10l [TTJ [12] . In 
SET regime a NEMS can be represented as a SET device coupled to a mechanical 
(harmonic) oscillator. The coupling is provided by a force (typically of electrostatic 
origin, see e.g. Ref. [13]) acting on the oscillator. The value of the force depends on the 
charge state of the SET device, providing the feedback. This feedback can be strong 
even if the electron-phonon coupling is weak: Indeed, if one consideres an SET next to 
the Coulomb blockade threshold, only two charge states, say, n = and n — 1, with n 
being the number of electrons at SET, are important. The occupation of SET above the 
Coulomb blockade threshold fluctuates between zero and one, providing a force which 
is a random telegraph signal. The signal swings an underdamped oscillator to large 
amplitudes, which in their turn affect the transport properties of the SET device by 
modifying the tunnel matrix elements and the energy difference between the two charge 
states. 

The transport properties of NEMS devices depend on a number of parameters. One 
is the coupling strength g, defined as the square of the ratio of the displacement of the 
oscillator center under the action of the force at the amplitude of the zero-point motion, 
g = F 2 /(hMu 3 ), where M and u are the mass and the frequency of the oscillator, and 
F is the difference of the forces acting on the oscillator in the charge states n — 1 and 
n = 0. In is important that the Planck constant is in the denominator, and thus strong 
coupling actually means quantum regime. Other parameters determining the behavior 
are given by ratios of various energy scales characterizing the NEMS device. The most 
relevant one is the ratio between the typical tunnel rate T and the oscillator frequency 

UJ. 

Surprisingly, the strong-coupling (quantum regime) is more extensively studied in 
the literature. Assuming the oscillator is underdamped, for g > 1 and u 3> T, the 
behavior of the system is dominated by Franck-Condon physics. The tunnel rates are 
modified due to emission of phonons in the course of electron tunneling [3]. For the 
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coupling to a single mode, this leads to the steps in the current as the function of 
the applied bias voltage, the height of the step is determined by the coupling constant 
[TU [15l [16]. For u <C T the electron level width becomes bigger than the distance 
between the steps, and Franck-Condon structure disappears. However, in this regime 
the motion of the oscillator is much slower than the electron tunneling, and one can 
use Born-Oppenheimer approximation, considering electrons in the quasi-stationary 
field potential provided by the oscillator. In this situation, the oscillator may become 
bistable, and the electron tunneling is dominated by switching events between the two 
states of the oscillator [T7J 13 [19] . All these phenomena are the manifestation of strong 
mechanical feedback. Similar effects have been studied in Ref. [23] in the context of 
superconducting NEMS. 

It is less obvious that strong feedback is also possible at weak coupling g <C 1. This 
is the classical regime, where Boltzmann equation serves as the starting point [20| 121] . 
Ref. [22] studied the low-frequency case w C T. The oscillator motion in this case 
is described by Fokker-Planck equation with effective diffusion and effective damping 
(originating from the electron tunneling out of the SET device), both determined by 
the energy dependence of the tunnel rates. It turns out that the strong feedback regime 
is feasible, but the behavior of the current strongly depends on the energy dependence 
of the tunnel rates. For instance, in the two most commonly investigated cases - 
electron tunneling through a single level and electron tunneling through a continuum of 
levels with the constant density, like in a single electron transistor, the strong feedback 
does not appear. Four distinct regimes have been identified: (i) no phonons generated; 
(ii) the oscillations are generated with the fixed finite amplitude; (iii) the oscillator 
is bistable, one state has the oscillations with zero amplitude, and another one the 
oscillations with a finite amplitude; (iv) the system is bistable, with the two states 
representing the oscillations with two different amplitudes. In the regimes (ii), (iii), 
and (iv) current is strongly modified with respect to the case when the phonons are 
not generated. A quantity even more sensitive to the oscillations is the current noise. 
The natural measure of the current noise in SET devices is Poisson value of the zero- 
frequency spectral density S, Sp = 2el, I being the average current [21]. It turns out 
that the noise is strongly enhanced avove the Poisson value in the regimes (ii), (iii), and 
(iv), and may even become super-Poissonian in the regime (i), when the current is not 
renormalized. 

In this Article, we consider the only regime not addressed so far: weak coupling 
g < 1 and high frequency uj ^> V. We show that it shares many features with the 
low-frequency classical regime described above. Even though we have not been able to 
identify any bistable regimes, we still find the two regimes of zero and finite amplitude, 
the latter demonstrating strong mechanical feedback. We find that the noise, in contrast 
to the low-frequency case, is always sub-Poissonian. 
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2. Boltzmann equation 

At weak coupling, the motion of the oscillator is classical, and the behavior of the 
system is characterized by the joint distribution function, P n (x, v,t). Here, n is the 
charge state of the SET. We assume that the SET is biased close to the edge of one of 
the Coulomb diamonds, so that only two charge states are important for transport, for 
definiteness n = and n = 1. Furthermore, x and v are the coordinate and the velocity 
of the oscillator. The starting point of our classical approach is the Boltzmann (master) 
equation for the distribution function [20, [22]. 



T = - M^x ~^ + F„; (2) 



St [P] = (2n - 1) (T+(z)P„ - r-(x)Pi) , (3) 

which holds for an arbitrary relation between uj and V . Here, the total force T acting 
on the oscillator is the sum of the elastic force, friction force, and charge-dependent 
coupling force, respective to the order of terms in Eq. Q 1 is the quality factor. 
We count the position of the oscillator from its equilibrium position in the n = state. 
In this case, F n = nF . The "collision integral" St [P] in the right-hand side represents 
single electron tunneling. There are four tunnel rates, R , where the subscripts L 
and R denote tunneling through the left or right junction, and the superscripts + and 
— correspond to the tunneling to and from the island, respectively; = + T R . 
Each rate is a function of the corresponding energy cost AE L R associated with the 
addition/removal of an electron to/from the island in the state n = 0/1 via left or right 
junction (AE^ R = —AEj jR ). Two independent energy differences are determined by 
electrostatics and depend linearly on the voltages. There is also a contribution to each 
energy from the shift of the oscillator, AE^ = AE^ ^ — Fx, AE R = AE R ^ +Fx, where 
AE° are the corresponding energy differences in the absence of mechanical motion. 

The condition uj ^> T t = T + +r~ means that the motion of the oscillator is very slow 
compared with the typical time an electron spends in the SET device. In this situation, 
the probabilities Pq and Pi average over the fast osicllator motion. In the leading order, 



if we parameterize x = esinyj, v = cue cos cp, e = y2E/Mu 2 , these probabilities do not 
depend on (p. Consequently, we expand the probabilities in the following way, 

P„(x,v,t) w P%(E,t) + cos <pF*(E,t) + sin <pPZ(E,t) , 

with P^' c P®. Here and below we disregard the terms proportional to smmip, cosmip, 
with m > 2. This procedure is similar to the transformation of Boltzmann equation 
into the diffusion equation in the semi-classical theory of electron transport in metals. 

We can obtain a closed set of equations for P® ,b,c by multiplying Eq. (TjQ) with 1, 
cos ip and siny? and subsequently averaging over the phase, throwing out sin2yj and 
cos2(/? terms. It is important that the tunnel rates only depend on the coordinate and 
not of velocity of the oscillator, and thus are functions of simp and not cos p. After 
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some algebra, we obtain Pq l = Pq = 0, 

PS(E,t) = <£lp(E,t), Pi(E,t) = J^P(E,t), 

(It) (sin y?r ) 

where the angluar brackets denote averaging over the phase (p, and the function P(E) 
obeys the equation 




ghu 3 E (r+> 



2 (r t ) 

_dP / 1 (siny?r+)(r-) - (siny?r-)(r+) 
\~dE + V 2ghuE (r+)(sin 2 ^T 3 ) 

and we have disregared the mechanical damping Q -1 - 



3. Distribution and current 

The stationary solution of Eq. (J4]) is easily found as 



P(E) = P(0) exp (- j(E')dE' 



(5) 



l{E) 



1 (sin <pT+) {T-} - (sin ^r~) (r+) 



I 2ghuE (r+)(sinVr-) 
Note that, similarly to the low- frequency case u <^T t , the distribution (jSJ) is very 
sharp. Indeed, the typical energy E is of the order of the applied voltage eV, and all 
energies which occur in Eq. (jSJ) are in our classical consideration much bigger that 
the quantum energy fkv of the oscillator. The function 7 can have positive as well as 
negative sign. If it becomes negative at some values of energy (amplification instead of 
the dissipation), there is a possibility that the strong feedback regime emerges. 

To illustrate the energy dependence of the probability, we have chosen exponential 
energy dependence used previously for low frequencies in Ref. [22] . 

It* = 2Vl R e- a ^ AE U(l - f F (-AE+ R )) ; 

rz,* = Tl R e a ^ AE ^ M AE IjR ) , (6) 

the factor 2 accounting for the spin degeneracy of the state n = 1. For concrete 
illustration, we choose = 0.3 , a« = 0.75 , and T° L = T° R . In the figures, W = eV g , the 
gate V g and bias V voltages are measured in the units of W c , the parameter associated 
with the energy dependence of the tunneling rates, with the value smaller than the 
charging energy. 

Fig. 1 presents the regions in gate-bias voltage plane corresponding to the two 
regimes. In the regime (i), the probability distribution has a sharp peak around 
zero energy. Since the energy represents the mechanical motion of the oscillator, it 
corresponds to the positive damping at any energy and the absence of mechanical 



Mechanical feedback in the high-frequency limit 



6 




-4-2 2 

W/W c 



Figure 1. The stability regions in the gate-bias voltage plane. Bold solid lines indicate 
the edge of the Coulomb diamonds. Insets show the sketch of P(E) in each region. 
The horizontal lines indicate bias voltages used for current and noise scans in Figs. 2 
and 3. 



feedback — phonons are not generated by the electron tunneling events. In the regime 
(ii), the distribution function sharply peaks around a finite value of the oscillator energy: 
Mechanical oscillations are generated. We were not able to detect the existence of the 
bistable regimes, similar to (iii) and (iv) described in the Introduction. The regime (ii) 
of phonon generation only emerges outside the Coulomb diamonds. 
For (i), the distribution function can be approximated as 

P(£)= 7 (0)exp(- 7 (0)£) , (7) 

where we have normalized the solution. For (ii), we have the Gaussian centered around 
the most probable value E m , 




The current is found as 

I(E)P{E)dE, 7(^7) = e <EiKS>Z^S><£l> . ( 9) 

Fig. 2 shows the results for the voltage dependence of the current. The trace (a) is taken 
in the regime (i), and the current is not modified by mechanical motion. The traces 
(b), (c) and (d) cross the region (ii), and the current dependence in this regime deviates 
from the one without mechanical motion. The deviations are stronger for higher bias 
voltages; additional peaks in the current develop. 

4. Current noise 

The current noise spectral power at zero frequency is found from the expression [22] 

/*oo /*oo 

S = -A SI(E)u(E)dE, 51(E) = 1(E) - / I(E)P(E)dE , (10) 
Jo Jo 
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Figure 2. Current modification in strong feedback regime for different bias voltages. 
The dashed (solid) lines give the current modified (unmodified) by mechanical motion. 
The modification is restricted to the regime to region (ii) where the phonon generation 
of takes place. 



and u solves the equation Cu = SI(E)P(E). In contrast to the low-frequency case, £ is a 
first-order differential operator. This fact simplifies the calculations and readily provides 
analytical estimates for noise in the two regimes. We parameterize u(E) = v(E)P(E), 
and the equation for v reads 



lghuj 3 E (T+) dv 



51(E) . 



V 2 (r t ) dE 

Solving it and substituting the result to the expression for noise, we obtain 



S 



ghixr Jo 



5I{E)P{E) j* ^5I(E')dE' . 



This expression can be evaluated with the use of the approximations flj]), flS}. 

In the regime (i), we write 51(E) « (gtuv /2e 2 )d 2 1 / 'dV 2 (E — 7 _1 (0)), where the 
second derivative is evaluated at E = 0. Substituting this into the expression for noise, 
we obtain 

S 



'-(g^f 1 



d 2 I 



e 2 dVf, 



<r t > 



E=0 



(T+> 



(7(0)) 



-5/2 



(12) 



The comparison with the Poisson value of noise Sp = 2el gives the estimate S/Sp ~ 
(T 2 /to 2 )(ghuj/eV g ) 3 ^ 2 . In the high-frequency regime T t <C to, both factors are small, and 
thus in the regime (i) the mechanically-induced noise is always sub-Poissonian. 
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(a)eV/W c =1.1 (b) ev/W c =2.2 
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Figure 3. Mechanical contribution to current noise for different bias voltages, (a) 
region (i); (b) region (ii). 



In the regime (ii), expanding SI 
16 \2EZ ( dl \ 



3a; V ghuj \ edV g 



{dI{E m )/dE){E - E m ), we obtain 
-l 



E—E„ 




E—E„ 



<r+> 



(13) 



The estimate for the noise power is S/Sp ~ (T 2 t / uo 2 )(ghoj / eV g ) l l 2 . It is not surprising 
that the noise is relatively higher than in the regime (i), where there is no mechanical 
motion induced. However, the noise is still below the Poisson value, which means that 
the behavior of the current noise is dominated by the shot noise. 



5. Conclusions 



In this Article, we studied current and current noise for an SET device coupled to 
an underdamped harmonic oscillator in the only regime not considered so far: weak 
coupling ^ C 1 and high frequency uj ^> Y t . We find that, similarly to other regimes, 
coupling to mechanical modes of the oscillator - phonons - excited by the tunneling of 
electrons through the SET device, may have a strong effect on the transport properties 
of the SET device. This is the phenomenon of strong mechanical feedback. 

However, we also find that the feedback effects are the weakest of all the regimes. 
It was very much expected that the effect of phonons on electron transport is stronger 
for stronger coupling. But we also find that for high frequencies the effects are 
less pronounced than for low frequencies. Indeed, no bistability regimes have been 
discovered, the strong feedback only manifests in oscillations with a fixed amplitude. 
The mechanical contribution to the current noise is small as compared with the shot 
noise contribution. It still can be serparated from the white shot noise due to its 
frequency dependence, however, the noise behavior is less spectacular than for low 
oscillator frequencies, where it is sometimes exponentially enhanced in comparison with 
the Poisson value. 

In this Article, we assume that the energy dependence of the tunnel rates is not 
too strong. All our results are expressed in terms of the tunnel rates averaged over the 
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period of mechanical oscillations. One can imagine an opposite situation — strongly 
energy-dependent rates (but at each energy still below the frequency of the oscillator). 
Since the oscillator moves quickly on the scale of the typical time the electron spends in 
one of the SET charge states, the instant tunnel rate performs fast oscillations with a big 
amplitude. In this situation, it is the easiest for an electron to tunnel when the tunnel 
rate is the highest, which typically would correspond to the maximum displacement of 
the oscillator. In this synchronization regime, the electron jumps are synchronized with 
the oscillator period. The condition for the appearance of the syncrtonization regime 
is (dT/dE)5E 3> T, where 5E ~ (eV) 2 / (ghuo) is the shift of the mechanical energy 
if the oscillator is displaced between the two extreme positions. Presumably, in the 
syncronization regime the transport properties are determined by the maximum (rather 
than the averaged) value of the tunnel rate over the oscillation period. Detailed analysis 
of the sincronization regime lies outside the scope of this Article, but we expect that the 
current is enhanced as compared with the "regular" strong-feedback regime considered 
above, whereas the current noise is suppressed since the electron stream becomes more 
regular. 
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